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Quantum coherence is a fundamental common trait of quantum phenomena, from the interference of matter 
waves to quantum degeneracy of identical particles. Despite its importance, estimating and measuring quan¬ 
tum coherence in generic, mixed many-body quantum states remains a formidable challenge, with fundamental 
implications in areas as broad as quantum condensed matter, quantum information, quantum metrology and 
quantum biology. Here we provide a quantitative definition of the variance of quantum coherent fluctuations 
(the quantum variance) of any observable on generic quantum states. The quantum variance generalizes the 
concept of thermal de Broglie wavelength (for the position of a free quantum particle) to the space of eigenval¬ 
ues of any observable, quantifying the degree of coherent delocalization in that space. The quantum variance is 
generically measurable and computable as the difference between the static fluctuations and the static suscepti¬ 
bility of the observable; despite its simplicity, it is found to provide a tight lower bound to most widely accepted 
estimators of “quantumness” of observables (both as a feature as well as a resource), such as the Wigner-Yanase 
skew information and the quantum Fisher information. When considering bipartite fluctuations in an extended 
quantum system, the quantum variance expresses genuine quantum correlations (of discord type) among the two 
parts. In the case of many-body systems it is found to obey an area law at finite temperature, extending therefore 
area laws of entanglement and quantum fluctuations of pure states to the mixed-state context. Hence the quan¬ 
tum variance paves the way to the measurement of macroscopic quantum coherence and quantum correlations 
in most complex quantum systems. 

PACS numbers: 03.65.Ud, 42.50.Lc, 05.30.-d 


I. INTRODUCTION 

Quantum mechanics introduces two fundamentally new 
traits in a physical system: 1) an intrinsic uncertainty on the 
knowledge of observables (Heisenberg’s uncertainty or coher¬ 
ent quantum fluctuations), and 2) a superior form of corre¬ 
lation among degrees of freedom, stemming from correlated 
quantum uncertainties (or entanglement) [1-3]. Quantum un¬ 
certainty of observables persists even at zero temperature in 
the form of so-called zero-point fluctuations, responsible for 
macroscopic quantum phenomena such as the inability of liq¬ 
uid Helium to solidify at ambient pressure [4]. On the other 
hand two quantum systems (hereafter called A and B) can 
exhibit correlations far exceeding any classical counterpart, 
which for pure quantum states are embodied by entanglement 
[5]. The supremacy of both fluctuations and correlations in 
quantum systems, as compared to classical ones, is at the heart 
of the complexity of many-body quantum states, challenging 
all realms of quantum physics, from relativistic quantum held 
theory to atomic/molecular physics and quantum condensed 
matter. At the same time quantum fluctuations and correla¬ 
tions (going beyond entanglement [6]) are by now recognized 
as essential ingredients for the supremacy of quantum devices 
over classical ones, in the context of both quantum informa¬ 
tion processing [6] and quantum metrology [7]. 

Despite their fundamental as well as practical importance, 
quantum coherence and quantum correlations remain very 
hard to both quantify theoretically and to measure experi¬ 
mentally. Quantum uncertainty of an observable and quan¬ 
tum entanglement between two subsystems are generically 
defined only for pure states [1, 5]. In the case of generic. 


real-life mixed states, the most widespread concept of quan¬ 
tum coherence is related to the thermal de Broglie wavelength 
(TdBWL) [8], expressing the spatial extent of coherent quan¬ 
tum fluctuations for a single quantum particle in free space; 
but this concept does not even extend to systems as simple 
as a particle in a potential. More recently, several definitions 
of mixed-state quantum coherence have been put forward [9- 
14], which nonetheless share the generically prohibitive re¬ 
quirement of knowing the full density matrix of the state, 
and they are therefore limited to few-body systems. As for 
the entanglement of mixed states, one can only provide suffi¬ 
cient conditions (witnesses) for the presence of entanglement 
between the components of the system [15-17]. Yet, even 
for non-entangled mixed states it has been recognized that 
quantum correlations may exist, associated with the violation 
of classical information-theory identities, and quantified via 
the so-called quantum discord [6] and discord-like quantities 
[11, 18]. Despite their deep conceptual meaning, entangle¬ 
ment and quantum correlations of mixed states remain in gen¬ 
eral information-theoretical objects, generically accessible (to 
calculations and measurements alike) only when defined be¬ 
tween two (or a few) elementary quantum degrees of freedom 
[18-23]. 


Here we show that both quantum coherent fluctuations and 
quantum correlations in generic quantum states can be quan¬ 
tified in terms of elementary physical concepts. The variance 
of fluctuations in generic mixed states possesses in fact an ad¬ 
ditive structure, in which the incoherent/thermal part can be 
separated from the coherent/quantum part, which we name 
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FIG. 1: Thermal vs. quantum fluctuations, (a) Different imaginary-time paths 0(r) in the space of eigenvalues of the observable O are 
shown, associated with the path-integral representation of a generic mixed state p. While the thermal/incoherent fluctuations {5'^0)t are 
associated with the fluctuations of the centroids of the paths (dashed blue lines), the quantum/coherent fluctuations {5^0)q are associated with 
the fluctuations of the paths around their centroids; (b) Geometry of the A-B bipartition of an extended quantum system used in the text; (c) 
Scaling of the total (tot), thermal (T) and quantum (Q) fluctuations of the staggered particle number Ns ,a on a subsystem A of size La for 
hardcore bosons on a square lattice at temperature T/J — 0.5 (the system is defined on a L x L torus with L = 32). All fluctuation terms 
exhibit volume-law scaling. Here and in the following graphs, the error bar is smaller than or comparable to the line thickness; (d) Temperature 
dependence of the Ns ,a fluctuations for a subsystem of linear size La = L/2. The dashed line indicates the infinite temperature limit, in 
which each lattice site fluctuates independently, with a shot-noise variance n(l — n) where n = 1/2 is the lattice filling. 


quantum variance (QV). The QV is defined in terms of the 
violation of a classical, static fluctuation-dissipation relation, 
and as such it is fully computable and measurable for generic 
systems. The QV of a given observable is a measure of its 
genuine quantum uncertainty in mixed states, and, in the case 
of bipartite fluctuations, it represents a measure of correlated 
quantum uncertainties, namely of quantum correlations. Re¬ 
markably the QV is convex (namely it decreases upon in¬ 
coherent mixing of states with the same QV), and it gives a 
tight lower bound to both the Wigner-Yanase skew informa¬ 
tion [9, 10] and to the quantum Fisher information [24] which 
are widely accepted, yet generically prohibitive measures (for 
both calculations and experiments) of the quantumness of ob¬ 
servables and of correlations [11, 12, 18, 25-27]. 


The structure of the paper is as follows. Sec. II intro¬ 
duces the separation between coherent and incoherent fluc¬ 
tuations, and the definition of QV; Sec. Ill reviews the fun¬ 
damental properties of the QV as a lower bound to known 
measures of “quantumness” of observables; Sec. IV describes 
the volume-law scaling of QV for generic extensive observ¬ 
ables; Sec. V illustrates the fundamental separation of scales 
between thermal and quantum fluctuations at a thermal criti¬ 
cal point; Sec. VI discusses the area-law scaling of the QV of 
bipartite fluctuations; Sec. VII illustrates the link between the 
QV and other measures of quantum correlations in the case 
of free fermions; and finally Sec. VIII elaborates on the gen¬ 
eral link between QV on the one side, and entanglement and 
quantum correlations on the other; and on the experimental 
measurement of the QV with cold-atom quantum simulators. 
The technical aspects are kept to a minimum level in the main 
text, and they are postponed to the Appendices. 


II. SEPARATING CLASSICAL AND QUANTUM 
LLUCTUATIONS. 

Let us first show that, given a density matrix p such that 
(...) = Tr[/3(...)]/Tr/5, and a generic Hermitian operator O, 
the fluctuations of the latter can be written as 

{5^6) = (6^) - {Of = {6^d)T + {s^d)Q ( 1 ) 

where {5‘^0)t represents thermal/incoherent fluctuations, 
while {5‘^0)q represent quantum/coherent fluctuations. In 
the following we shall focus our attention on thermal equilib¬ 
rium states, but the whole treatment is readily generalizable 
to arbitrary density matrices (see App. B). If p = e~^'^/Z 
(Z — Tr(e“^^)) is the thermal density matrix of a system of 
Hamiltonian TL at temperature ksT = 1/P, and = 0, 

it is well known that the fluctuations of O satisfy a (classical) 
fluctuation-dissipation theorem 

(5^0) = xo keT ( 2 ) 

where xo = d^F/dh^\h=o is the susceptibility associated 
with the application of a term —hO to the Hamiltonian, and 
F — —ksTlogZ is the free energy. On the other hand, if 
\0, TL] ^ Q the quantum uncertainty on the value of O adds 
up to the thermal fluctuations, and, as a result 

{5‘^0) > xo ksT =]- f dr (jd(T)jd(O)) =: 

P Jo 

( 3 ) 

where 0(r) = is the operator evolved in imagi¬ 

nary time. Eq. (3) shows that thermal fluctuations do not ex¬ 
haust the total fluctuations of the observable. It is then natural 
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to define the QV for the observable O as the residual fluctua¬ 
tions, or as the violation of the classical fluctuation-dissipation 
relation of Eq. (2): 

{6^6)Q = {6^6)-xokBT. (4) 

The QV has a particularly suggestive interpretation in the con¬ 
text of a path-integral representation of the partition function 
of the system, using a basis of Hilbert space which diagonal¬ 
izes the O operator (see Fig. 1(a)). As discussed in App. B, 
this allows one to cast the partition function in the form: 

where 0{t) is a periodic trajectory (0(0) = 0(^)) in the 
space of eigenvalues of O, associated with the Feynman path 
in the basis diagonalizing O, and S is the associated action 
weighting the trajectory. When assigning a path-integral ex¬ 
pression to each of the terms in Eq. (4), one can easily find 
that (see App. B) 

{S^d)Q = (^^JdT{0{r)-0f'^^ ( 6 ) 

where (...)s is the average over the ensemble of paths 0(t) 
weighed by the action S, and 

O = 0[0{t)] = ^ J dr 0{t) (7) 

is the centroid of the path [28]. Eq. (6) shows that the QV rep¬ 
resents the (squared) amplitude of the imaginary-time fluctua¬ 
tions of the trajectory 0(t) around the path centroid. Clearly 
such fluctuations have a genuine quantum origin [28-30]. If 
O is the position i of a one-dimensional particle, in App. C 
we show that {6'^x)q ^ A^g, namely the QV is tightly related 
to the quantum uncertainty on the position expressed by the 
TdBWL AdB- When moving to higher dimensions and generic 
quantum systems, the QV generalizes therefore the concept of 
TdBWL (or quantum coherence length) to the space of eigen¬ 
values of any Hermitian operator. And, most remarkably, it 
does so in a computable and measurable manner, being ex¬ 
pressed as the difference between a fluctuation property and a 
response function. 


III. PROPERTIES OE THE QUANTUM VARIANCE. 

The QV represents a physically measurable lower bound to 
fundamental quantities in quantum information. The Dyson- 
Wigner-Yanase skew information [31] 

4(0,p) = -^Tr{[d,p“][d,pi-“]} (8) 

with a G [0,1], generalizing the Wigner-Yanase skew infor¬ 
mation (a = 1/2) [9], probes the quantum uncertainty of O 


stemming from its non-commutativity with p. As shown in 
App. D, the QV is simply 

(< 5 ^ 0 ) q [ p ] = [ da la- ( 9 ) 

Jo 

From the convexity of la [31] follows the convexity of the 
QV. Moreover one can prove that 

(<520)q[p]</i/2(6,p) (10) 

(the equality holding for pure states). Finally, the quantum 
Fisher information [24] 

Fq{0',P) = 2^ |(n|0|TO)|^(Pn -Pm)^/(Pn +Pm) (H) 

nm 

(where p = Pn/n) (n|) expresses the sensitivity of the den¬ 
sity matrix to a unitary transformation U (h) = gener¬ 

ated by the observable O, and it quantifies the fundamental 
metrological gain in using the state p to estimate the parame¬ 
ter h [7]. As shown in App. E, 

{S^d)Q[p]<FQ{d;p)/d. (12) 

The inequality becomes an equality for pure states. As dis¬ 

cussed later, the inequalities satisfied by the QV have consid¬ 
erable implications concerning its importance for entangle¬ 
ment witnessing and metrological applications. Conversely, 
the computability and measurability of QV gives unprece¬ 
dented insight into the skew and quantum Fisher information 
for quantum many-body systems. 

IV. QUANTUM VARIANCE OE A GUOBAE OBSERVABLE 
AND VOLUME LAW. 

Due to its inherent quantum nature, the QV exhibits very 
special size and temperature dependences. In the following 
we shall concentrate on thermal equilibrium states, and we 
start our analysis with the case of a generic, macroscopic ob¬ 
servable O that does not commute with the Hamiltonian of 
the system. As an example we consider the case of two- 
dimensional hardcore bosons on the square lattice: 

•H =-J^ (SJSj+h.c.) (13) 

(ij) 

where bi, S| are hardcore boson operators, satisfying the 
(anti)-commutation relations {bijbj} = 1 and [bi,bj] = 
[bi, 6^] = 0 (i ^ j). We treat this model with a numerically 
exact quantum Monte Carlo algorithm based on the Stochastic 
Series Expansion approach [32], which allows us to investi¬ 
gate the imaginary-time dynamics of many-body observables 
[33]. The Hamiltonian Ti does not commute with any finite- 
wavevector Fourier component of the density profile, and in 
particular with the staggered particle number 

Ns = J2i-iybik- 


(14) 
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FIG. 2: Critical thermal fluctuations, non-critical quantum fluc¬ 
tuations. (a) Temperature dependence of the order parameter fluc¬ 
tuations at the Ising transition of 2d hardcore bosons with nearest 
neighbor repulsion V = 2.1J; the sharp peak marks the transition at 
Tc ~ 0.78 J; (b) Scaling of fluctuations close to the critical tempera¬ 
ture: the total and thermal fluctuations are found to scale as 
with 7 = 7/4 and v = 1 for the 2d Ising universality class. 


To investigate the scaling of fluctuations (both thermal and 
quantum) we isolate a subsystem A of linear size La in a 
larger system (of linear size L - see Fig. 1(b)), and we inves¬ 
tigate the scaling of local observables/fluctuations in A with 
the size of the A region itself. This approach allows one to 
extract scaling properties while using a single simulation box, 
and it is also directly applicable to experiments giving access 
to local properties, such as those based on quantum-gas mi¬ 
croscopy [34]. 

Fig. 1 (c) shows that both the thermal and the quantum con¬ 
tribution to fluctuations obey a volume-law scaling in the ex¬ 
ample at hand: 

{5^Ns,a)t, {5^Ns,a)q ~ ( 15 ) 

where iVg a = ^ volume-law scaling of 

quantum fluctuations is generically expected when the observ¬ 
able of interest is extensive, and its Heisenberg’s uncertainty 
is the result of the non-commutativity between an extensive 
set of terms in the Hamiltonian and in the observable in ques¬ 
tion. The separation between thermal and quantum fluctua¬ 
tions gives rise to a very meaningful result when tracking the 
temperature dependence of the fluctuations on a subsystem of 
fixed size (La = T/2 in this case). As shown in Fig. 1(c), the 
thermal component grows linearly with T at low T, while the 
quantum component decreases monotonically with T starting 
from the zero-point fluctuations. Most remarkably, in the ex¬ 
ample at hand quantum fluctuations are found to dominate the 
total fluctuations, and they lead to a monotonically decreasing 
behavior of {6‘^Ns), in complete contradiction with the classi¬ 
cal expectation that fluctuations should grow with temperature 
at low T. 


V. QUANTUM VARIANCE DOES NOT GO CRITICAL AT 
A THERMAL TRANSITION. 

Having shown that QV generically obeys a volume law for 
extensive non-conserved observables, one can naturally ask 
what is the fate of QV at a thermal critical point, at which 
thermal fluctuations of the order parameter become super- 
extensive. If the QV only captures the quantum mechanical 
part of fluctuations of the order parameter, one would natu¬ 
rally expect that its scaling is not modified at a thermal tran¬ 
sition, given that the latter is purely driven by thermal fluc¬ 
tuations. To answer to this question, we consider a quantum 
many-body model exhibiting a thermal phase transition with 
an order parameter not commuting with the Hamiltonian; this 
is readily obtained by generalizing the hardcore-boson Hamil¬ 
tonian of Eq. (13) to include a nearest-neighbor repulsion V: 

H = -JJ2 (blbj + h.c.) (fij - . 

(U> (U> ^ ^ ^ ^ 

( 16 ) 

When V > 2J the model has an Ising phase transition at fi¬ 
nite temperature, marking the onset of a checkerboard density 
wave, with an order parameter given by the staggered density 
Ng. Hence, as in the previous section, it is meaningful to in¬ 
vestigate the temperature and size scaling of the fluctuations 
of the local staggered density Ng. In particular, to mimic the 
behavior in the thermodynamic limit (in which (iVg a) 0), 
we focus on the fluctuations around a finite-size estimate of 
the order parameter in the symmetry-breaking (SB) phase, 
given by (|iVs,A|): 

(S^Ng^A)^^^') = (NIa) - {\Ns,A\f ■ ( 17 ) 

Fig. 2(a) shows that the total and thermal fluctuations of the 
order parameter exhibit a sharp peak at the Ising transition 
temperature, while the QV is very smooth at the transition. In 
particular, the QV is a monotonically decreasing function of 
temperature, dramatically showing that quantum fluctuations 
cannot grow under incoherent thermal mixing of the density 
matrix, even when such fluctuations become singular at a crit¬ 
ical point. The monotonicity with temperature is in fact a 
different condition - required by physical arguments - than 
the convexity discussed in Sec III, which only accounts for 
the evolution under linear mixing. A closeup on the scaling 
close to the critical point (Fig. 2(b)) finds that the total and 
thermal fluctuations exhibit the critical super-extensive scal¬ 
ing , where 7 and v are the critical exponents for the 

susceptibility and correlation length. On the other hand the 
quantum variance maintains a volume-law scaling as in the 
non-critical regime. Therefore a critical point marks a net sep¬ 
aration of scales between thermal and quantum fluctuations 
of the order parameter, the latter being essentially irrelevant 
in the thermodynamic limit. This observation substantiates 
the common wisdom that quantum mechanics is irrelevant for 
the universal properties at thermal critical points, and it shows 
that order parameters close to a critical point have the nature 
of emergent classical observables. 
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FIG. 3: Bipartite fluctuations, (a)-(c) Scaling of local particle-number fluctuations in a subsystem A for square-lattice hardcore bosons 
{V = 0) at three different temperatures (T/J = 0.01, 0.1 and 0.4). Other parameters as in Fig. 1(c); (d) Temperature dependence of the 
particle-number fluctuations; same parameters as in Fig. 1(d). 


VI. QUANTUM VARIANCE OF BIPARTITE 
FLUCTUATIONS AND AREA LAW. 

The scaling properties of the QV change drastically when 
considering the case of bipartite fluctuations of an otherwise 
globally conserved quantity, such as the particle number TV. 
Such fluctuations have been the subject of several recent stud¬ 
ies in view of their relationship to entanglement in the case 
of pure states [35-38], as well as at finite temperature, for 
which a suggestion of how to extract the quantum contri¬ 
bution to fluctuations has been proposed in Ref. [36]. In 
the case of mixed states, [N, TL] = 0 implies automatically 
that {6'^N)q = 0. Taking then any bipartition of the sys¬ 
tem into A and B subsystems, imaginary-time fluctuations 
of the local particle numbers Na and TVs are perfectly an¬ 
ticorrelated, so that the QV in each subsystem is the same, 
{S'^Na)q = {6'^Nb)q- Perfect correlation in the quantum 
uncertainties of Na and Nb suggests that the QV captures 
genuine quantum correlations between A and B whenever ap¬ 
plied to bipartite fluctuations of globally conserved quantities. 
Remarkably Fig. 3(a-c) shows that the QV of bipartite fluc¬ 
tuations scales like the boundary between A and B, thereby 
obeying a so-called area law 

{S^Na)q - L‘^a~" , (18) 

namely the extensive (volume-law) part of bipartite fluctua¬ 
tions is entirely of incoherent origin. This strongly suggests 
that the QV captures the fluctuations associated with coher¬ 
ent particle exchanges at the boundary between A and B. For 
the hardcore-boson problem at hand, such fluctuations obey a 
logarithmically corrected area law at T = 0 (when all fluctu¬ 
ations are quantum) [37, 38], 

{S^Na)q^ L<^A~"^ogLA (19) 

turning then into an area-law scaling at finite T. Nonethe¬ 
less a logarithmic violation can still be observed at sufficiently 
low temperature and for small sizes of the A region - namely 
smaller than the thermal correlation length ^ for density fluc¬ 
tuations [56]. Interestingly, the area-law scaling of the QV 



FIG. 4: Quantum correlations vs. quantum mutual information. 

(a) Scaling of the quantum variance of bipartite particle-number fluc¬ 
tuations and of the quantum mutual information in a system of free 
fermions ona L x L square lattice (L — 32) at half filling for three 
different temperatures (J is the hopping amplitude); (b) temperature 
dependence of the same two quantities, along with the total entropy 
Sa, the total fluctuations {S'^Na), the Wigner-Yanase skew informa¬ 
tion I\/ 2 {Na, p), and the quantum Fisher information Fq{Na', p), 
for an A region with linear size La = L/2. The T~^ decay of the 
mutual information at high temperature has been proven rigorously 
for free fermions in Ref. [39], and it is proven for the quantum vari¬ 
ance, skew information and quantum Fisher information in App. F. 


(either straight or logarithmically violated) is found to domi¬ 
nate the scaling of total fluctuations at sufficiently small sizes 
La of the subsystem A, as shown in Fig. 3(a-c). This makes 
the (logarithmically violated) area law of bipartite quantum 
fluctuations observable under experimentally realistic condi¬ 
tions. 
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VII. QUANTUM VARIANCE PROVIDES QUANTUM 
CORRELATIONS. 

The area-law scaling of bipartite QV of mixed states sug¬ 
gests a link to the similar scaling exhibited by entanglement 
entropy in ground states of local Hamiltonians [40]. Yet 
measures of entanglement at finite temperature (such as the 
negativity [16]) do not admit a simple physical interpreta¬ 
tion in terms of entropy of quantum fluctuations (but see be¬ 
low for further discussion on QV and entanglement). As al¬ 
ready pointed out, QV is rather connected to quantum corre¬ 
lations, a more general concept than entanglement. The ex¬ 
istence of such correlations is captured by the quantum dis¬ 
cord [6], given by the difference between the quantum mu¬ 
tual information, I{A : B) = Sa + Sb — Sab (or the en¬ 
tropy sub-extensivity due to correlations between A and B) 
and the classical mutual information J{A : B) (or the maxi¬ 
mum amount of information gained on A by performing mea¬ 
surements on B). Here Sa(b) = ~^'^Pa{b) ^ogpA{B) is the 
entropy of the reduced density matrix of subsystem A{B), and 
Sab is the total entropy. The operation of maximization im¬ 
plicit in the definition of quantum discord makes it generically 
non-computable when A and B are extended subsystems of a 
quantum many-body system. 

On the other hand, in some special systems the existence of 
quantum correlations is witnessed by more accessible quan¬ 
tities. Indeed we argue that, in the case of an ideal gas, any 
form of correlation stems from the quantum statistics, while 
it is trivially absent in the classical limit. The existence of 
correlations between A and B is generically captured by the 
quantum mutual information, whose nonzero value is then a 
direct proof of quantum correlations existing in the system. 
[57] In the case of an ideal lattice gas the existence of corre¬ 
lations between A and B stems physically from the coherent 
exchange of particles at the A-B boundary, and hence it is 
tightly linked to the quantum fluctuations of particle numbers. 
In the following we shall particularly focus on the case of a 
free Fermi gas on a lattice at half filling, for which the mu¬ 
tual information and QV of particle-number fluctuations can 
be easily calculated via exact diagonalization [41]. 

The quantum mutual information of free fermions ex¬ 
hibits an area law at finite temperature, as recently proven in 
Ref [39, 42]. Fig. 4 shows that the area law of mutual infor¬ 
mation and of QV of particle-number fluctuations are related, 
as the prefactors of the thermal area laws, I {A : B)/2 « 

ai{T)LA and {6‘^Na)q « aAi{T)LA, are proportional at 
2 

all T, a I « ^uat. Remarkably, this is the same relation¬ 
ship holding between the particle-number variance and the 
entanglement entropy for free fermions at T = 0 (Fig. 4(a)) 
[36, 43], and between total entropy and variance in a degener¬ 
ate Fermi gas (Fig. 4(b)). Hence the particle-number QV pro¬ 
vides experimental access to the mutual information of free 
fermions at finite T, as much as the total variance of particle- 
number fluctuations gives access to the entanglement entropy 
in the ground state. Moreover, as shown in Fig. 4(b) the QV 
provides a meaningful lower bound to both the Wigner-Yanase 
skew information and to the quantum Fisher information; 


in particular at high temperatures we And that « 

lli/ 2 iNA,p) and {S^Na)q « |(FQ(iV^; p)/4) . 

VIII. CONCLUSIONS AND OUTLOOK. 

In conclusion, we have introduced the quantum variance 
of generic observables, generalizing the concept of quantum 
Heisenberg uncertainty to the case of mixed states - and acting 
as the “thermal de Broglie wavelength” in the space of eigen¬ 
values of arbitrary observables. In the case of bipartite fluctu¬ 
ations, the QV expresses the quantum correlations among the 
two subsystems in arbitrary mixed states. The quantum uncer¬ 
tainty may dominate the fluctuations in quantum many-body 
systems, leading to a completely non-classical behavior (fluc¬ 
tuations decreasing with temperature, scaling of fluctuations 
obeying area laws or logarithmically violated area laws, etc.). 

Owing to its definition in terms of fully measurable quanti¬ 
ties (fluctuations and response function, see Eq. (4)), the QV 
is readily accessible to state-of-the-art experiments. All the 
requirements for the measurement of the QV, and in particu¬ 
lar of its scaling in a bipartite setting, are met by trapped-ion 
experiments [44] as well as quantum-gas microscope experi¬ 
ments [34], enabling access to local degrees of freedom. As 
an example, in microscopy experiments recent progress [45] 
has demonstrated the ability to resolve different single-site oc¬ 
cupation numbers {n = 0,.., 3) in an optical lattice, providing 
access to local fluctuations. Moreover the local response func¬ 
tion can be probed by an increase of the local uniform (or stag¬ 
gered) chemical potential in a given region A of the system, 
making use of holographic masks [34] (see App. G for further 
details). Hence the total and thermal fluctuations are indepen¬ 
dently accessible, as well as the scaling of their difference (the 
QV) with subsystem size (Figs, l(c-d), 3 and 4). 

At the theory level, the QV represents a most accessible 
way to assess quantum correlations, entanglement, and the 
metrological use of quantum many-body states. As proposed 
in Refs. [11, 18], the minimal skew information and quan¬ 
tum Fisher information associated with local observables in a 
subsystem A are both discord-type measures of quantum cor¬ 
relations, and the latter dictates the minimal precision on the 
estimation of the parameter of an arbitrary local unitary oper¬ 
ation; the QV offers a natural measurable lower bound to both 
quantities (see App. D and E for further discussion). More¬ 
over both the skew information [25] and the quantum Fisher 
information [26, 27] of collective spin variables witness en¬ 
tanglement among k qubits when exceeding a fc-dependent 
bound; a similar violation of the bound by the QV is therefore 
an even stronger witness - see App. H for a detailed discus¬ 
sion. 

The QV lends itself to analytical as well as to large-scale 
numerical simulations based e.g. on quantum Monte Carlo - 
as shown in the present work. Hence its study can be readily 
extended to generic quantum many-body systems at equilib¬ 
rium, including interacting fermions, quantum spin models, 
etc. as well as to non-equilibrium mixed states. While we 
have mostly focused our attention on bipartite correlations, 
an extension of our study to multipartite correlations can also 
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be readily achieved by introducing the concept of quantum 
covariance, as we will further develop in future work. This 
opens the perspective of a complete characterization of quan¬ 
tum correlations in extended quantum systems, based on ex¬ 
perimentally accessible quantities. 
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Appendix A: Operator approach to coherent vs. incoherent 

fluctuations in mixed states: area law of quantum coherence 

As pointed out in the main text, in the case of mixed states 
described by a density matrix p there is a fundamental dis¬ 
tinction between thermal/incoherent fluctuations and quan¬ 
tum/coherent fluctuations of any observable which does not 
commute with p. This distinction is best captured via the 
path-integral representation of the density matrix, as discussed 
in the following App. B. Here we give an alternative picture 
solely based on the operator picture of the density matrix. In 
the following we shall choose, as observable of interest, the 
particle number Na in the region A of the system, capturing 
the quantum correlations between the region in question and 
its complement. 

When leaving the ground state of local Hamiltonians, one 
expects to encounter states with the generic feature of possess¬ 
ing volume-law entanglement, and volume-law fluctuations of 
particle number [46]. Hence one may naively suspect that, 
when dealing with excited states, the quantum coherent fluc¬ 
tuations are stronger, and not weaker, than in the ground state. 
This is indeed true, but it is only meaningful provided that, 
in an experiment, one is able to deterministically prepare one 
and the same excited state, in order to probe its fluctuation 
properties over many shots of the experiment itself. This last 
requirement is generally prohibitive, as experiments on ex¬ 
cited states generally probe the properties of ensembles (every 
shot of the experiment reproducing a different state). Whence 
the relevance of the concept of density matrices p, not only in 
the context of systems coupled to dissipative baths, but also 
in the context of systems evolving uniquely under their own 
Hamiltonian dynamics. 

In the latter case, let |'I'(f)) be the instantaneous state of the 
system, and let 0 be a time window sufficiently long for time 
averages to equal ensemble averages (namely averages over 
repeated shots of the experiment). Then the density matrix 
describing the ensemble is well described by 

p^^j%t\^{t)){^{t)\ (Al) 


Despite the fact that each state may exhibit volume-law 

entanglement and coherent fluctuations, the ensemble proper¬ 
ties are quite different. Indeed, we can write |'I'(f)) as 

(A2) 

Na {rai}iv^ 

where |{ni}; Na) is aFockstate {rii} characterized by having 
Na particle in A. After time/ensemble averaging, the density 
matrix takes the form p = pb + Pod where 

PD = (A3) 

E E E Na) {{n'i}, Na\ 

Na WAnj^ 

is the diagonal part of the density matrix (in terms of the quan¬ 
tum number Na), and 


POD = (A4) 

E E E P{ni},NA',{n'A,N'j^ I ^ a) ({^i}) ^ A I 

Na^N'^ {ni}^^ {^[liv^ 


is the off-diagonal part; here 


1 / 

P{ni},NA-,{n'J,N'^ = Q I (f) . 

(A5) 

It is evident that [pu,^a] = 0 while [pod,-A^] 4^ 0. 
Therefore the off-diagonal part, containing the coherence be¬ 
tween configurations differing by the number of particles in 
A, is the part of p responsible for the quantum fluctuations 
of Na, captured by the quantum variance. The total, ex¬ 
tensive fluctuations of Na are given by the diagonal part, 
{5^Na) = {Nl)p, - {Na)% ; as a consequence the quan¬ 
tum coherent contribution, which stems from the off-diagonal 
terms, remains hidden in this calculation, and it cannot be for¬ 
mally separated from the incoherent part. The proper separa¬ 
tion between incoherent and coherent fluctuations is achieved 
within the path-integral formalism, as described in the main 
text and below in App. B. Nonetheless the operator form of 
the density matrix provides further insight into the physical 
origin and spatial structure of coherent quantum fluctuations, 
as discussed below. 

The instantaneous coherence (t) con¬ 
nects states with Na — N'j^ ^ as it is typical of ex¬ 

cited states in Hilbert space. But the time/ensemble-averaged 
coherence P{ni},NA-,{n'.},N'j^ in Efl- (A5) has a much shorter 
range away from the diagonal. Indeed, assuming that {ui} 
and {n'} are connected by moving m particles from sites 
Jij ■■■, jm to sites ii,..., im, one has 

P{ni},NA-,{n'J,N'^ = (A6) 


Tr 




\{ni},NA){{ni},NA\ 


where bi, bl are the destruction/creation operators of the par¬ 
ticles of interest (of arbitrary statistics). Hence, as one may 
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have expected, the magnitude of P{ni},NA-,W},N'^ is con- where 0{t) is the continuum limit of the sequence 
trolled by that of the 2TO-point correlation function, namely {Oai, Oa^-, OomI’ 


\P{n,},NA-,{n',},N'^ \ < l(^]i ’ ’ ' ■ (A7) 

In general such a correlation function will exhibit a fast decay 
with the (minimum) distances between pairs of points ip and 
jq. This in turn implies that, in order to have a sizable co¬ 
herence (Eq. (A5)), two configurations {ui} and {n'} should 
differ by particle moves which, when occurring between A 
and its complement B, are localized (algebraically or expo¬ 
nentially) around the boundary between the two regions. This 
observation generally excludes a volume law for the coher¬ 
ent part of particle-number fluctuations, and it leaves an area 
law (up to multiplicative logarithmic corrections) as the only 
possibility. 


Appendix B: Path-integral representation of a generic density 
matrix and of the quantum variance 

In this section we derive the path-integral representation for 
a generic density matrix, generalizing the discussion of the 
main text to arbitrary mixed states beyond thermal equilib¬ 
rium. Moreover we derive the path-integral expression for the 
quantum variance. 

Any (semi-positive definite) density matrix p can always be 
cast in the form 


M-l 


■,-S 


M—¥oo ' ^ 




(B4) 

Once the density matrix has been given the thermal form 
Eq. (Bl), it is straightforward to deform the density matrix 
upon application of a held h coupling to O, 


p{h) 


^-PCH-hd) 

Tr[e-/3(«-'iO)] ’ 


(B5) 


which allows one to define the response function in the stan¬ 
dard way as xo = JrTr[/5(/i)6] 

h—0 

The path-integral representation of response function leads 
to the expression 


Xo = 


dr 50 {t)50(Q)^ 

— (^J dr J dr' 50 {t)50(t') 


(B6) 


where 50{t) = 0{t) — {0)s, and we have invoked the peri¬ 
odicity of 0{t) paths in imaginary time. Here 


{■■■)s = ^[ V[0{t)] (...) e-^ (B7) 

^ JO{0)=O{p) 


e-/3« 

Tr[e-/3«] ’ 


(Bl) 


namely in the form of a thermal density matrix with (effec¬ 
tive) temperature ksT — 1//3. Eor generic (non-thermal) 
mixed states the specific value of (3 is completely irrelevant, 
and one could set in the following /3 = 1 in some convenient 
energy units; yet, in order to make contact with the case of 
thermal equilibrium, hereafter we will keep the inverse tem¬ 
perature /? explicitly indicated. We consider a generic observ¬ 
able O which is diagonalized by a basis |Oq, {A}q,), where 
Oa is the eigenvalue for O, and {A}q, are the other quantum 
numbers possibly labeling the state. The partition function 
Z = Ti[exp{—(did)] can be cast in the form of the trace of the 
product of infinitesimal propagators between successive states 
|Oai,{A}aJ, namely 


M-l 

M—^oo 

{a.} 

(B2) 

where X]{a } ^ short-hand notation for the multiple sum 

over the quantum numbers (0„. ,{A}q,J labeling each state 
in the propagation sequence ai,a 2 , ...jaM = oti- Summing 
over the A quantum numbers, and taking the limit M —>■ oo, 
one obtains formally the path-integral expression 

Z = f V[0(t)] e-^^Oir),d^o{T),...] (23) 

Jo{o)=oiP) 


is the average over the space of 0(t) paths. Moreover one has 
that 


(S^O) = ^QdT . (B8) 


Combing Eqs. (B6) and (B8), one readily obtains the path- 
integral expression for the quantum variance 


{S^O)q 


((5^0) - xo ksT 



0{t) 


1 


2 


dr' 0 (t') 


(B9) 

s 


showing that it represents the average variance of fluctuations 
of 0(t) paths around their centroid. 

We end this section by noticing that the deformation of the 
density matrix to Eq. B5 is a physically meaningful opera¬ 
tion for thermal states - as it can be obtained by turning on 
the perturbation —hO in the Hamiltonian within an isother¬ 
mal setting - see App. G for further discussion in the specific 
case of quantum gas microscopes. Hence in the case of ther¬ 
mal states, neither the measurement nor the calculation of the 
quantum variance requires the full knowledge of the density 
matrix. On the other hand, for generic mixed states the defor¬ 
mation of p should be thought of in general as a mathematical 
operation. Devising physical (namely, experimentally realis¬ 
tic) operations that can lead to the deformation of a generic 
density matrix as in Eq. B5 is an outstanding task, which we 
postpone to future investigations. 
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60 


function of temperature. In the case of the free particle the frequency 
LU is to be understood as an arbitrary constant setting the energy scale 
hu) and length scale Uho- 


Appendix C: Quantum variance as generalized de Broglie 
wavelength 

1. Quantum variance for simple models 

The analytical calculation of the quantum variance of the 
position operator is illuminating in the case of simple mod¬ 
els, namely the free particle and the harmonic oscillator in 
one dimension. It is convenient to start from the second one, 
and to obtain the free-particle result as a limiting case. In the 
case of the harmonic oscillator, the position fluctuations are 
readily obtained from the diagonal part of the density matrix 
{x\e~^'^\x)/Z [28], while the susceptibility Xx = to a 
displacement of the harmonic oscillator potential ^muix'^ 
\mujx^ — hx is readily obtained by the linear displacement of 
the average, {x) —>■ (x) — h/imuP'). As a result the quantum 
variance takes the form 


FIG. 6: Quantum variance as Heisenberg’s uncertainty for mixed 
states, (a) The path-integral representation of the partition function 
for a single one-dimensional particle at T —^ 0 (or for a generic 
pure state) is a sum over infinitely long paths in imaginary time 
x(t); on average each path visits the region lx,x + dx] a number 
of times proportional to \'tj){x)\'^dx, where ij){x) is the ground-state 
(or, more generally, the pure-state) wavefunction. The width of the 
wave function’s square modulus gives the Heisenberg’s uncertainty, 
namely the amplitude of coherent quantum fluctuations; (b) In the 
case of mixed states of a generic many-body system, the probability 
distribution for the coherent quantum fluctuations of a generic ob¬ 
servable O is instead given by the probability that the (imaginary- 
time) instantaneous value of 0(r) differs from the path centroid 
O = f dTO{T). 


to the case of a particle in an external potential, such as the 
case of the harmonic oscillator. In Fig. 5 the quantum vari¬ 
ance of the position for the two models discussed above shows 
the expected monotonic decrease with temperature, due to the 
shrinking of the imaginary-time “duration” of coherent quan¬ 
tum fluctuations. 




ajo \ sinh(l/6l) 

2 _cosh(l/61) — 1 


(harm, osc.) 

(Cl) 


where Oho = \/hf(jruJ) and 9 = kBT/huj. In the limit T —i- 
0 one recovers Heisenberg’s uncertainty in the ground state. 


{5‘^x)q = alJ2. 

On the other hand, taking the limit w —?> 0 gives the re¬ 
sult for the free particle, which, after careful expansion of 
Eq. (Cl), gives 


{5‘^x)q = ^ XIb{T) (free particle) (C2) 


where XdBiT) = /{rnksT) is the thermal de Broglie 

wavelength. The link between the quantum variance and the 
de Broglie wavelength shows that the quantum variance of the 
particle position gives the characteristic (squared) amplitude 
of coherent quantum fluctuations at finite temperature [28]. 
In fact one may interpret the quantum variance of the posi¬ 
tion as a generalization of the thermal de Broglie wavelength 


2. Imaginary-time fluctuations as quantum coherent 
fluctuations 

As seen in App. B, the quantum variance of a generic ob¬ 
servable gives the characteristic amplitude of fluctuations for 
such an observable along the imaginary-time dynamics of the 
system. On the other hand, in App. C we have established a 
direct relationship between the quantum variance of a free par¬ 
ticle and the thermal de Broglie wavelength, namely the char¬ 
acteristic width of wavepackets at finite temperature. In this 
section we bring the two observations together to argue that 
the quantum variance generalizes the concept of thermal de 
Broglie wavelength, or finite-temperature coherence length, 
to the space of eigenvalues of any observable (not only the 
position operator) and for any quantum system. 

The relationship between the quantum variance and the 
de Broglie wavelength is very natural when considering the 
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fundamental link existing between the wavefunction of a 
pure state and the statistics of Feynman paths. Thanks to 
the parametrization in Eq. (Bl), a pure state can always be 
thought of as the /3 —?> oo limit of a mixed-state density ma¬ 
trix, and hence represented in the form of a path integral. In 
the case of a one-dimensional quantum particle, the path in¬ 
tegral for a pure state with wavefunction runs over in¬ 
finitely long trajectories x{t), whose fluctuations 6x around 
the average position (x) = {ilj\x\tjj) obey the statistics dictated 
by the modulus of the square function [47], namely, taking an 
arbitrary time r: 

P{Sx) ~ ^ J <5(a;(T) — (x) — Sx) e~^ 

= li’iix) + Sx)\^ . (C3) 

In particular any infinite trajectory contributing to the path in¬ 
tegral has the same statistical properties as the whole ensem¬ 
ble, so that the centroid of the path x = a;[a:(T)] must cor¬ 
respond to the expectation value (x). Hence, as depicted in 
Fig. 6(a), the imaginary-time fluctuations span the support of 
the wavefunction, and the (squared) amplitude of fluctuations 
of the path x{t) around its centroid x - the quantum variance 
- is the same as the (squared) width of the wavefunction, giv¬ 
ing Heisenberg’s uncertainty. In the case of free particles, the 
thermal de Broglie wavelength generalizes Heisenberg’s un¬ 
certainty on the position to the case of thermal states. There¬ 
fore, it is not too surprising that the quantum variance follows 
the de Broglie wavelength at finite temperature, as shown in 
App. C. 

The concept of quantum variance extends all the above con¬ 
siderations to generic observables and generic quantum sys¬ 
tems. The quantum variance provides the width of the prob¬ 
ability distribution for the fluctuations of generic observables 
around the path centroid (see Fig. 6(b)), namely 

p{SO) = ^ J 'D[0{t)] S{0{t) -0-60) (C4) 

where 0 = 0[0{t)] is the centroid. As seen in the case of 
the position of a one-dimensional particle, for a pure state the 
width of P{60) expresses the Heisenberg’s uncertainty on the 
observable O. When applied to a mixed state, the quantum 
variance generalizes therefore Heisenberg’s uncertainty, ex¬ 
pressing the (squared) amplitude of coherent quantum fluctu¬ 
ations of the observable. 


Appendix D: Quantum variance vs. skew information and local 
quantum uncertainty 

In this section we shall discuss the relationship between the 
quantum variance and the skew information [9], the latter be¬ 
ing a widespread concept in quantum information to quantify 
the quantum uncertainty of an observable. In particular we 
shall show that the quantum variance provides a tight lower 
bound, based on physical observables, to the otherwise ab¬ 
stract skew information. Moreover the discussion of the rela¬ 
tionship between the skew information and the quantum vari¬ 
ance allows one to conclude on the convexity of the latter. 


Finally we will see how the quantum variance relates to the 
recently introduced “local quantum uncertainty” [11], which 
is advocated as a measure of quantum correlations. 


1. Wigner, Dyson, Lieb and the convexity of quantum variance 

The Dyson-Wigner-Yanase (DWY) skew information [9, 

31] 

4(0,/3) = -^Tr{[6,p“][d,pi-“]} (0<a<l) (Dl) 

probes the quantum uncertainty of the observable O due to 
its non-commutativity with the density matrix of the system. 
Replacing O by 60 = O — (O) does not alter the above defi¬ 
nition. 

Writing again the generic density matrix p as a thermal 
state, p = e~^'^IZ (with an arbitrary effective inverse tem¬ 
perature /3), one can immediately show that the DWY skew 
information can be expressed as an imaginary-time correla¬ 
tion function 

/„(0, p) = (6^0) - (<56(r = a(3)66{Q)) . (D2) 

Hence, clearly, the DWY skew information 1^ expresses 
the amount by which the imaginary-time correlation function 
{60{t)60{0)) at a time r = aP has decreased with respect to 
the equal-time (t = 0) value. Hence the DWY skew informa¬ 
tion probes the imaginary time fluctuations in a similar man¬ 
ner to quantum variance. As a consequence a link between 
the two quantities can be expected, and it is straightforwardly 
established in the form 

{6^d)Q = {6^0) -^fdr {6dir)6dm 

da 4(0, p) (D3) 

namely the quantum variance is equal to the average DWY 
skew information. In particular the Wigner-Yanase (WY) 
skew information [9], Ia^ij 2 , is an upper bound to the DWY 
skew information 



4(6,p)<4^i/2(6,p) (D4) 

as it is easy to prove due to the monotonic decay of imaginary- 
time correlation functions, {60{t)60{0)) > {60{t = 

(6/2)60(0)). As a consequence one readily obtains that the 
quantum variance is always lower than the WY skew infor¬ 
mation 


(<J'0)Q<4=i/2(d,p). (D5) 

Lieb [48] proved that the DWY skew information is convex 
for any value of a, namely: 


4(0, Aipi -f A 2 P 2 ) < Ai4(0, pi) + A24(0, P 2 ) (D6) 
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for any real numbers Ai, A 2 . Using Eq. (D3), the property of 
convexity is immediately inherited by the quantum variance. 
The convexity of quantum variance is a fundamental figure of 
merit to assess the quantum variance as a probe of quantum 
coherent fluctuations: if pi and p 2 have the same quantum 
variance, any linear incoherent superposition of the two has 
necessarily a lower quantum variance. 

2. Quantum variance vs. local quantnm uncertainty 

Given an A-B bipartition of an extended quantum system, 
Ref [11] has introduced the concept of local quantum uncer¬ 
tainty (LQU) 

ip) = minQA / 1/2 (Oa, p) (D7) 

as the minimum of the WY skew information over all 
local observables 0\ in A having a given spectrum 
A. Ref. [11] argues that this observable-independent (but 
spectrum-dependent) quantity acts as a discord-like quantity, 
namely a measure of quantum correlations between A and B. 
In order to capture quantum correlations among all, equally 
weighted degrees of freedom of A and those of B, it is ob¬ 
vious to request that the observable Oa be an extensive one, 
namely the sum of local observables Oj, 

^^macro) = ^ . (D8) 

j&A 

This ensures that quantum fluctuations of all degrees of free¬ 
dom in A are taken into account into the skew informa¬ 
tion. As the size of A grows, the spectrum of the operator 
»^icro) _ j becomes a continuous one, and it 

is contained in a finite interval [Amin > Amax] ■ This behavior 
applies to all extensive operators of the kind of ^ ^nd 

their spectrum can easily be reduced to one and the same A by 
a simple shift and rescaling in the definition of the operator. 
Hence, in the sense of Ref. [11], one can define a macroscopic 
LQU defined as a minimum over all operators 

O^acro)^ which is arguably the most appropriate definition 
of discord-type correlations among all degrees of freedom of 
A and those in B. We assume that A and B interact with a 
coherent Hamiltonian term, leading to an exchange of energy, 
and possibly also particle, or magnetization, etc... Hence in 
the minimization procedure we explicitly exclude the possi¬ 
ble existence of local conserved quantities ^ 

which would trivially lead to a vanishing macroscopic LQU. 

It follows immediately from Eq. (D5) that the macroscopic 
LQU is lower-bounded by the minimum quantum variance of 
macroscopic observables 

^^macro)(^) ^ min^™ / 1 / 2 p) 

> min^^rnacro) )q . (D9) 

The minimization implied by Eq. (D9) is readily performed 
for the quantum variance: the minimum quantum variance 


of macroscopic observables is realized by bipartite fluctua¬ 
tions of an otherwise conserved quantity, namely Oa such that 
\Oa + Ob: p] =0. Lor general quantum systems, the above 
requirement applies to the local energy, and, in the presence of 
a continuous symmetry, to the local particle number (for par¬ 
ticle models) or to the local magnetization (for spin models), 
etc. - assuming that the latter quantities are not conserved. 
In the case of equilibrium states of local Hamiltonians, we 
have shown in this work that the quantum variance of bipar¬ 
tite fluctuations obeys an area law: as a consequence, Eq. (D9) 
implies that the macroscopic LQU obeys at least an area law. 
On the other hand, in the ground state the WY skew informa¬ 
tion reduces to the variance of the operator 

h/2[0\ )> (DIO) 

and the scaling of the minimum variance of local macroscopic 
operators in the ground state of local Hamiltonians satisfies a 
(logarithmically violated) area law [38]. Even though the 
temperature dependence of the WY skew information is not 
generally known in the literature, one can assume that it is 
maximized at T = 0 (this is the case of the free-fermion ex¬ 
ample studied explicitly in Sec VII). Under this assumption, 
and given that the WY skew information at T = 0 coincides 
with the total variance, we obtain the inequalities 

min^(macro) ((5^)Q (Dll) 

< < min^^„acro) (J2o(macro)^(y ^ q) 

implying that the macroscopic LQU obeys at most a logarith¬ 
mically violated area law, namely 

\ogLA) . (D12) 

Appendix E: Quantum variance vs. quantum Fisher 
information: quantum correlations and metrology 

In this section we focus on the relationship between the 
quantum variance and the quantum Lisher information [24], 
a central quantity in quantum metrology due to its link with 
the maximum precision achievable in the estimation of the 
parameter of a given unitary transformation. Similarly to the 
skew information, the quantum variance offers a lower bound 
to the quantum Lisher information; we shall exploit this fact in 
the context of the recently introduced “interferometric power” 
[18] to explore the importance of the quantum variance of bi¬ 
partite fluctuations both for metrology and for quantum cor¬ 
relations. Lurther implications of this bound in the context of 
entanglement witnessing will be discussed in Sec H. 

1. Quantum variance as a lower bound to the quantum Fisher 
information 

The quantum Fisher information (QFI) [24] expresses the 
“distinguishability” (in the sense of the Bures distance) be¬ 
tween two density matrices p(h) and p{h -f 5h), belonging 
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to a family p{h) continuously parametrized by the parameter 
h. If the family of density matrices is obtained via a uni¬ 
tary transformation generated by an Hermitian operator O, 
p(h) = e~''^^p{h = 0)6**^^, the QFI takes the explicit form 

Fq{0,p) = EnmG'F(p„,Pm) \{n\5d\m)\^ where and 
\n) are eigenvalues and eigenvectors of the density matrix, and 


(El) 

Pn + Pm 

This is to be compared with the expression 
of the quantum variance, namely {S^O)q = 
Enm V {Pn.Pni)\{n\SO\m)\^ whcrc 


^ Q V {Pn : Pm ) 


Pn Pm 
2 


Pn Pm 

log(p«) - log(Pm) 


(E2) 


Comparing the two functions it is easy to realize that 

> Gqy^x, y) Q<x,y<\ (E3) 


whence the announced inequality 


concept of interferometric power (IP) of an observable in 
a bipartite {A + B) system with density matrix p as 

^Aip) = ]^^^Fq{Oa,p) (E8) 

4 oi 

where Fq{0\^p) is the quantum Eisher information associ¬ 
ated with a unitary transformation generated by a local ob¬ 
servable 0\ acting on A, and with spectrum A. The IP has 
a direct metrological meaning: it expresses the worst-case- 
scenario uncertainty (in the sense of the Cramer-Rao bound 
[50]) that one can achieve in the estimation of the parameter 
of a unitary transformation generated by an arbitrary observ¬ 
able which is local in A and has a given spectrum A. Ref. [18] 
argues that the IP is another discord-type measure of quantum 
correlations between A and B, leading to the conclusion that 
quantum correlations are a resource for metrology. 

It is immediate to see that the above conclusions carry au¬ 
tomatically over to the case of the quantum variance. Using 
the inequality Eq. (E4), one immediately has that 

Viip) > mm{S^d^)Q . (E9) 

oi 


> {S^d)Q . (E4) 

Incidentally we notice that the WY skew infor¬ 
mation admits a similar expression Ii/ 2 { 0 ,p) = 

E„m Gi^^^{pn,Pm) |(n|50|m)^ with 

Gl^/^{Pn,Pm) = - y/PnPm ■ (E5) 

Direct inspection into the G functions reveals the inequality 
chain: 

> Gi^^^ {x, y) > Gqy{x, y) 0 < x, y < 1 

(E6) 

whence the ensuing hierarchy: 

> /i/2(6, p) > {d^d)Q . (E7) 

In particular we notice that QEI and WY skew information 
have often been invoked as sharing similar properties [49], 
yet a well-defined relationship as the inequality of Eq. (E7) 
between the two has not yet been established (in this respect. 
Ref [49] uses a QEI-like quantity, which is not equivalent to 
the actual definition of QEI). The quantum variance further 
explicits this relationship by finding a common, non trivial 
lower bound for both quantities. 


In App. D 2 we argued that macroscopic observables 
in A, having an extensive spectrum A, capture the quantum 
correlations between all degrees of freedom in A and those 
in B. In the case of such observables one can perform the 
minimization immediately for the right-hand side of Eq. (E9), 
identifying the operator with the one satisfying the 

condition _|_ ^^acro) ^ _ q (again, as in App. D 2 

we are excluding local conserved quantities from the mini¬ 
mization). Hence the quantum variance of bipartite fluctua¬ 
tions provides a lower bound on the IP of macroscopic ob¬ 
servables, and on the quantum correlations and metrological 
resource that this quantity expresses. Similarly to Eq. (D12), 
this lower bound allows one to establish an area law scal¬ 
ing (with at most logarithmic corrections) to the IP of macro¬ 
scopic observables, under the assumption (verified e.g. by free 
fermions as in Sec. VII) that the QEI is maximised at T = 0. 
In particular this bound is very instructive in terms of the 
metrological utility of many-body states: the maximum quan¬ 
tum variance of bipartite fluctuations, and hence the maximum 
IP, is achieved for states exhibiting power-law correlations, 
and specifically in the vicinity of quantum critical points - see 
also [51] for a recent calculation of the quantum Eisher infor¬ 
mation in exactly solvable models of quantum-critical points, 
which confirms this conclusion. 


Appendix F: Quantum variance, skew information and 
quantum Fisher information of bipartite fluctuations for free 
fermions 


2. Metrological implications: the interferometric power 


1. Quantum variance 


In a similar manner to the definition of local quantum un¬ 
certainty discussed in App. D2, Ref. [18] has introduced the 


In this section we calculate the quantum variance of local 
particle-number fluctuations in the case of free fermions on a 
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d-dimensional hypercubic lattice at half filling. The density- 
density correlation function is given by 

((5ni(T)(5nj(0)) = (FI) 

^ ^ (1 - fy) 

h,h' 

where fk = [exp(/?efe) -I- 1]“^ is the Fermi distribution, and 
Efc = —2 J J2a=x y cos{ka) is the dispersion relation. Inte¬ 
grating the correlation function to get (S'^Na) and {S‘^Na)t, 
one obtains the quantum variance in the form 

{S^Na)q = (F2) 

k.k' i,jGA 


which, when expanded at high temperature, leads to the be¬ 
havior 

Ii/2{Na,p)= (F7) 

k,k' iJ^A 

The above expression is very similar to Eq. (F4) for the quan¬ 
tum variance, confirming that the two quantities have the 
same high-temperature behavior (as well as the same zero- 
temperature value). In particular, when /I —0: 

{S^Na)q = I h/2{NA, P) + 0{I3^) . (F8) 

3. Quantum Fisher information 


tmEE'’"'"'"'’A( i-A-) 


In the high-temperature limit /? —> 0 the quantum variance 
reduces to 

{S^Na)q = (F3) 

§^E E .) (e,-e.0^+O(/33). 

k,k' i,jGA 


Finally, when considering the QFI for free fermions [51], 
one finds the following expression for bipartite particle- 
number fluctuations 


Fq{Na,p) = 

E E fk (1 - /fcO tanh^ 

i,jGA 


(F9) 


d(efc-efe/) 

2 


One observes that the term linear in /3 vanishes, so that the 
dominant temperature dependence goes like T“^. 

Assuming for A the geometry of a hypercube of side La, 
the double sum over the A region can be performed exactly, 
leading to 


which leads to the high-temperature behavior 

Fq{Na,p)= (FIO) 

k.k' i,j^A 


{6'^Na)q = 

+ 0{I3^) . 


(F4) 



sinc^[(/cQ, —fc^)LA/2] 
sinc^[{/cQ; —fc^)/2] 


) (efc - tk'Y 


Comparing again with Eq. (E4), one can conclude that, at high 
temperatures 


{5^Na)q = ^ + 0(/3") . (Ell) 


In the limit La —> oo one has that 


4. Discussion 


^ smc^[(fca - fca)L^/2] 
^ sinc^[(fca - fc[,)/2] 


OO 

2 5 {ka — k'^ — 27rn) . 

n——oo 


(F5) 

One sees immediately that this limit would highlight a term 
scaling as L\ (volume law), but that this term is actually van¬ 
ishing because = efc-i- 27 rne„ ■ Hence one is left with an 
area-law scaling term. 


2. Wigner-Yanase skew information 


Using Eq. (D2) leads immediately to the following expres¬ 
sion for the skew information of bipartite particle-number 
fluctuations of free fermions 


Ii/2{Na,p) — 


1 


k.k' i.jGA 


(F6) 


Hence, as anticipated in the main text, the quantum fluctu¬ 
ations captured by the quantum variance, the WY skew infor¬ 
mation or the QFI display the same high-temperature behav¬ 
ior up to a global prefactor. This leads to a coherent picture 
for bipartite quantum fluctuations of free fermions. While the 
calculation of the quantum variance is easily extended to arbi¬ 
trary many-body systems which can be treated with state-of- 
the-art numerics, the same is generally not true for the WY 
skew information nor the QFI - although, unlike the QFI, 
the WY skew information lends itself to path-integral Monte 
Carlo approaches probing imaginary-time correlation func¬ 
tions. On the experimental side, the WY skew information, 
being an imaginary-time correlation function, is not accessible 
to experiments as such. As for the QFI, Ref [51] has recently 
shown that it is potentially accessible to experiments when 
cast as a frequency integral involving the dynamic suscepti¬ 
bility; in this respect, the quantum variance has the advantage 
of being expressed solely in terms of static correlations and 
response functions. 
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and B, all incoherent fluctuations (either stemming from the 
grand-canonical ensemble or from other sources) are system¬ 
atically subtracted away in the quantum variance, if one is able 
to realize experimentally the deformation of the density ma¬ 
trix as in (B5). We argue that this is indeed the case when the 
total particle number obeys an arbitrary statistics, namely the 
case in which the density matrix takes the general form 


FIG. 7: Quantum-gas microscope setup to measure the quantum 
variance. Here we sketch a possible scheme to measure the quantum 
variance of the local particle number in region A by adding a box-like 
potential to a two-dimensional optical lattice. This potential induces 
a local increase in the chemical potential, allowing one to probe the 
response function as the particle-number difference between region 
Al and A\ supplementing this measurement with the one of particle- 
number variance in region A gives access to the quantum variance. 


Appendix G: Measurement of bipartite quantum variance with 
quantum-gas microscopes 


A concrete proposal to measure the quantum variance of 
bipartite particle-number fluctuations in the context of ultra¬ 
cold quantum gases is illustrated in Fig. 7. Recent progress 
in quantum gas microscopes allows one to trap atoms in sin¬ 
gle layers of an optical lattice with superimposed traps of ar¬ 
bitrary geometries [34] and to measure single-site particle- 
number occupations [45]. As already discussed in App. B, 
to access the quantum variance of the particle number Na in 
the subsystem A one needs to measure the total variance of 
fluctuations (S'^Na), as well as the response function 


XNa 


+ 5^a) ~ {NA){fJ-A) 
SfiA 


(Gl) 


where ^a is the local chemical potential in region A, cou¬ 
pling to the particle number Na- The two quantities {6'^Na) 
and xna need to be measured in the same conditions of tem¬ 
perature and (offset) chemical potential. A way to achieve 
this in cold-atom experiments is to use a “multiplexing” setup 
as in Fig. 7, in which one single trap geometry allows one 
to measure both quantities at once. Indeed monitoring fluc¬ 
tuations of Na in region A allows one to extract {6'^Na) the 
total variance; on the other hand, a box-like potential super¬ 
imposed to the optical lattice creates a local increase in the 
chemical potential, giving access to the response function as 
{{Na') — {NA))/SfiA- If the regions A and A! are built sym¬ 
metrically around the (global) trap center, and if thermal equi¬ 
librium is established across the system, one is ensured that 
the two quantities are measured in the same thermodynamic 
conditions of temperature and offset chemical potential. 

One may worry that in cold-atom experiments the total par¬ 
ticle number has wide shot-to-shot fluctuations going well be¬ 
yond a grand-canonical description, and that this may alter 
the estimate of the quantum variance, adding spurious contri¬ 
butions coming from experimental systematics. On the other 
hand, as discussed in the main text and in sections App. A 


P=\ ^Pexp(iV) V^e-^'^rN (G2) 

N 

where pexp{N) is the experimental particle-number statistics, 
accounting for systematic shot-to-shot fluctuations, and Vn is 
the projector onto the Fock subspace with N particles. The 
deformation of the Hamiltonian implied in Fig. 7 leads to the 
desired deformation of the density matrix; hence the quan¬ 
tum variance (and its peculiar size and temperature scaling) 
can be experimentally measured even without postselection of 
the measurement shots according to the total particle number, 
with the obvious caveat that one is not measuring properties 
of the grand-canonical ensemble but the ones of the artificial 
ensemble realized experimentally. 

A similar setup, and similar considerations, can be applied 
to measure the quantum variance of the staggered particle 
number. In that case, one needs to shine a weak superlattice 
potential with twice the lattice spacing of the primary poten¬ 
tial over the region A'. 


Appendix H: Quantum variance as multiparticle entanglement 
witness 

Let us consider a system of N qubits, with collective spin 
operators J = ^ I*® 

producible [25, 52, 53] if it can be written as 

\^k-p.od) = (HI) 

where is an (entangled) state of a block of Ni < k spins, 
with the constraint that Ni = N. A mixed state is then 
said to be /c-producible if it is an incoherent superposition of 
fcs-producible states with kg < k 

Pfc—prod — E Ps\'4^ks—prod) {'4^ks — prod I . (H2) 

S 

Using Eq. (E4) and the results of Refs. [26, 27], one can prove 
that for fc-producible states the quantum variance of the col¬ 
lective spin components j“, and the QFI associated to trans¬ 
formation generated by the J“, satisfy the inequality: 

(<5V“)q < FQ(J“;pfe_p,od)/4 < nk^ + {N - nkf (H3) 

where n = [N/k] is the integer part of N/k. In fact the ex¬ 
act same bound as in the last inequality of Eq. (H3) holds 
for the WY skew information [25], and again it carries over 
to the quantum variance thanks to the inequality in Eq. (D5). 
The inequality of Eq. (H3) can be readily generalized to more 
general degrees of freedom than qubits, namely to collective 
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operators C = where Ci is an operator with a bounded 

spectrum contained in the interval [cmin, Cmax] • In that case 
the inequality takes the form [7] 

4(^^C')q < Fq{C; pk-prod) 

< (Cmax - Cinin)^[nfc^ + (IV - nfc)^] . (H4) 

Hence, similarly to what was already found for the WY 
skew information [25] and the QFI [7, 26, 27], a violation of 
the inequalities in Eqs. (H3) or (H4) or for the quantum vari¬ 


ance is a strong indication of the existence of multiparticle 
entanglement among a least {k -f 1) degrees of freedom. The 
condition of violation is actually very strong, as the bound of 
the last inequality in is rather loose for thermal states. Indeed 
the bound of Eq. (H3) is valid for fc-producible pure states and 
mixed states alike, but, given that all the quantities in question 
(WY skew information, QFI, and quantum variance) are ex¬ 
pected to decrease under thermal mixing, the bound is much 
looser for thermal states, and all the more so the higher the 
temperature. 
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